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Method and System for Scheduling Cars in Elevator Systems 
Considering Existing and Future Passengers 

Field of the Invention 

[01] This invention relates generally to scheduling elevator cars, and 
more particularly to elevator scheduling methods that consider future 
passengers. 

Background of the Invention 

[02] Scheduling elevators in a large building is a well-known hard 
industrial problem. The problem is characterized by very large state 
spaces and significant uncertainty, see Barney, "Elevator Traffic 
Handbook,'' Spon Press, London, 2003. Typically, a passenger requests 
elevator service by pressing a call button. This causes the elevator 
scheduler to assign an elevator car to service the passenger. 

[03] The earliest elevator schedulers used the principle of collective 
group control. In this heuristic, the nearest car, in its current direction of 
travel, is assigned to service the passenger, see Strakosch, "Vertical 
transportation: elevators and escalators,'' John Wiley & Sons, Inc., 
1998. Such scheduling is sub-optimal and unpredictable. For this reason, 
collective control is unacceptable when passengers expect to be notified 
about which car will pick them up, immediately after the call is made. 
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[04] Another heuristic minimizes a remaining response time (RRT) for 
each passenger. The RRT defines the time it takes to pick up each 
passenger as prescribed by the current schedule, see U.S. Patent No., 
5,146,053, ''Elevator dispatching based on remaining response timej' 
issued to Powell et al., on September 8 1992. That heuristic focuses only 
on minimization the waiting time of passengers, and ignores altogether 
the effect of the current assignment on the waiting times of future 
passengers. 

[05] Within RRT-based minimization, a further distinction can be made 
between those methods that ignore the uncertainty associated with the 
desired destination floors of passengers, see Bao, ''Elevator dispatchers 
for down-peak traffic,'' Technical Report, University of Massachusetts, 
Department of Electrical and Computer Engineering, Amherst, 
Massachusetts, 1994, and those that properly determine the expected 
RRT of each passenger with respect to destinations, see Nikovski et al., 
"Decision-theoretic group elevator scheduling,'' 13* International 
Conference on Automated Planning and Scheduling, Trento, Italy, June 
2003, and U.S.Patent Application Sn. 10/161,304 "Method and System 
for Dynamic Programming of Elevators for Optimal Group Elevator 
Control," filed by Brand et al. on June 3, 2002, incorporated herein by 
reference. 

[06] However, the uncertainty associated with future passengers is 
entirely new matter for at least two reasons. Accounting properly for the 
effect of the current decision on the waiting times of all future passengers 
is an extremely complicated problem. First, the uncertainty associated 
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with future passengers is much higher because the arrival time, the arrival 
floor, and the destination floor are all unknown. Second, the current 
decision potentially influences the waiting times of passengers arbitrarily 
far into the future, which makes the theoretical optimization horizon of 
the problem infinite. 

[07] In spite of the computational difficulties, ignoring future 
passengers often leads to sub-optimal scheduling results. The current 
assignment affects the future movement of the cars, and influences their 
ability to serve future calls in the minimal amount of time. 

[08] One particular situation that exemplifies the importance of future 
passengers is peak traffic. During down-peak traffic periods, for example, 
at or near the end of the workday, most future passengers select the main 
floor as their destination. Because these future passengers are most likely 
distributed over upper floors, scheduling for down-peak traffic is a very 
hard problem. 

[09] During up-peak traffic periods, most future passengers arrive at the 
main floor and request service to upper floors. Typically, the up-peak 
period is much shorter, busier and concentrated than the down-peak 
period. Therefore, up-peak throughput is usually the limiting factor that 
determines whether an elevator system is adequate for a building. 
Therefore, optimizing the scheduling process for up-peak traffic is 
important. 
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[010] Consider the following scenario. A call is made at some upper 
floor. A single car is parked at the main floor, and the scheduler decides 
to serve the call with that car, based only on the projected waiting times 
of passengers. If the car at the main floor car is dispatched to serve the 
call, the main floor remains uncovered and future passengers will have to 
wait much longer than if the car had stayed. This shortsighted decision, 
commonly seen in conventional schedulers has an especially severe 
impact during up-peak traffic, because the main floor quickly fills with 
many waiting passengers, while the car services the lone passenger 
above. 

[Oil] Several elevator scheduling methods are known for considering 
future passengers, with varying success. Some schedulers use fiizzy rules 
to identify situations similar to the one discussed above and make 
decisions that are more sensitive to future events, see Ujihara et al., "The 
revolutionary AI-2000 elevator group-control system and the new 
intelligent option series,'' Mitsubishi Electric Advance, 45:5-8, 1988. 
However, that method has major disadvantages. First, the rules need to be 
coded manually. Therefore, the system is only as good as the 'expert'. 
Second the interpretation of fuzzy-rule inferences between the rules often 
behaves erratically, particularly when there is no applicable rule for some 
specific situation. Thus, the elevators often operate in an unintended and 
erratic manner. 

[012] Another method recognizes that group elevator scheduling is a 
sequential decision making problem. That method uses the Q-leaming 
algorithm to asynchronously update all future states of the elevator 
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system, see Crites et al., ''Elevator group control using multiple 
reinforcement learning agents,'' Machine Learning, 33:235, 1998. They 
dealt with the huge state space of the system by means of a neural 
network, which approximated the costs of all future states. Their 
approach shows significant promise. However, its computational 
demands render it completely impractical for commercial systems. It 
takes about 60,000 hours of simulated elevator operation for the method 
to converge for a single traffic profile, and the resulting reduction of 
waiting time with respect to other much faster algorithms was only 
2.65%, which does not justify its computational costs. 

[013] The prior art methods are either labor-intensive or computationally 
expensive or both. Therefore, there is a need for a method that optimally 
schedules elevator cars, while taking future passengers into 
consideration, particularly for up-peak traffic intervals. 

Summary of the Invention 

Brief Description of the Drawings 

[014] Figure 1 is a block diagram of an elevator system that uses the 
invention; 

[015] Figure 2 is a flow diagram of a method for scheduling elevator cars 
according to the invention; and 
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[016] Figure 3 is a grid showing Markov chains according to the 
invention. 

Detailed Description of the Preferred Embodiment 
System Structure 

[017] Figure 1 shows an elevator scheduler 200 according to our 
invention for a building 101 with upper floors 102, a main floor 103, 
elevator shafts 104, elevator cars 105. The main floor is often the ground 
or lobby floor, in other words the floor where most passengers entering 
the building mainly arrive. 

[018] For the purpose of our invention, passengers are formally classified 
into several classes according to variables that describe what is known 
about the passengers. The variables introduce uncertainty into the 
decision-making process of the elevator scheduler. The classes are riding, 
waiting, new and future passengers. 

[019] For each riding passenger 1 1 1, the arrival time, the arrival floor, 
and the destination floor are all known. The riding passengers are in cars, 
and no longer waiting. 

[020] For each waiting passengers 112, the arrival time, the arrival floor, 
and the direction of travel are known. The destination floor is not 
known. A car has been assigned to service each waiting passenger. 
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[021] For a new passenger 1 13, the arrival time, the arrival floor, and the 
direction of travel are known because the new passenger has signaled 120 
a call. The general problem is to assign a car to service the call of the new 
passenger. At any one time, there is only one new passenger. 

[022] The above three classes of passengers 111-113 are collectively 
existing passengers. The reason we call these passenger existing is 
because they have already arrived physically, and the system knows 
something about all of these passengers. Of the existing passengers, only 
the waiting passengers and the new passenger have non-zero waiting 
times. 

[023] For fiiture passengers 1 14, who do not exist yet, nothing is known. 
At best, the passenger variables can be described stochastically by 
random variables, or be estimated from past data. All passengers include 
existing and future passengers. 

[024] The specific problem is to assign a car to service the new 
passenger so that the expected waiting time for all passengers, existing 
and future, is minimized. 

Method of Operation 

[025] Figure 2 shows a method for scheduling cars of the elevator system 
100 according to the invention. The method 100 executes in response to 
a call 201. The call can be any floor. First, the scheduler 200 determines, 
for each car, based on future states 209 of the elevator system, a first 
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expected waiting time 21 1 for all existing passengers 1 1 1-1 13 if the car is 
assigned to service the call. Second, the scheduler determines, for each 
car based on a landing pattem 219 of the cars 105, a second expected 
waiting time 221 of the future passengers 1 14 if the car is assigned to 
service the call 102. For each car, the first and second expected waiting 
times are combined 230 to produce an adjusted waiting time 231, and the 
car with the lowest adjusted waiting time is assigned 240 to service the 
call 201. 

[026] Ideally, the elevator scheduler would determine the marginal costs 
of all possible assignments, with all sources of uncertainty integrated out, 
before making an assignment. However, due to the insurmountable 
computational complexity of the scheduling problem, the vast majority of 
commercial elevator schedulers typically resort to heuristic methods that 
ignore some or all of this uncertainty. 

[027] In typical up-peak traffic periods, a substantial number of future 
passengers, e.g., between 80% and 95%, arrive at the main floor. The 
waiting times of these main floor arrivals is the dominant component in 
the overall waiting time of an elevator system during up-peak traffic 
periods, and the current decision of an elevator scheduler should attempt 
to minimize the expected waiting time of passengers at the main floor. 

[028] Hence, we begin with a simplifying assumption that all future 
passengers arrive at the main floor. The effect of not modeling future 
arrivals at other floors shortens the time-horizon in which predicted waits 
are accurate to the near future. However, this effect is explicitly worked 
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into the calculations later, as a discounting factor. In addition, for up- 
peak traffic periods, most future passengers do in fact arrive at the main 
floor. 

[029] With this assumption, the current decision of the elevator scheduler 
affects the waiting times of future passengers through the future arrival of 
cars at the main floor. We call this sequence of arrivals of the cars at the 
main floor the landing pattern. 

[030] For the purpose of the invention, the landing pattem 219 of cars at 
the main floor is determined by the following factors. First, riding 
passengers at upper floors can select the main floor as their destination. 
Second, empty cars can automatically select the main floor as the place to 
park while waiting for a next call. Determining the landing pattem 219 
effectively marginalizes out individual future passengers 214. 

[031] Optimal parking strategies and their effects on the landing pattem 
are described in U.S. Patent Application Sn. 10/293,520 ''Optimal 
Parking of Free Cars in Elevator Group ControU^' filed by Brand et al., 
on November 13, 2002, incorporated herein by reference. 

[032] One strategy to service main floor passengers preferentially is to 
send each car to the main floor immediately after it has completed 
servicing the last riding passenger. For a building with C cars, a landing 
pattem 219 is an array of times 

T = [r;,r2, ...,rd, for 7,^0, 
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where 7} is the arrival time of car 7 = 1 , . . . , C at the main floor after it has 
delivered all of its riding passengers. 

[033] Because there is uncertainty about the destinations of the waiting 
passengers 1 12 and the new passenger 113, the landing pattem T is a 
vector- valued random variable with a probability distribution P(T), T ^ 
rover the space of all possible landing pattems T219. 

[034] Ideally, the scheduler 200 should determine an expected waiting 
time V(T) for each possible landing pattem T, and take the 
expectation of that time with respect to the probability distribution P(T) 

as ( P(T) >= L P(T)y(T)i/T. 

[035] Here ( > denotes the expectation operator. Indeed, this is an exact 
estimate of the waiting times of main floor passengers, under the above 
assumption that all new passenger arrivals are at the main floor. 
However, there is no practical way to determine the probability 
distribution P(T). Even if there was, the size of the space of all possible 
landing pattems is huge. Integrating over this space is computationally 
impractical. 

[036] Instead, we use a substitute landing pattem including individual 
expected arrival times at the main floor of each car T = [ ^ i, 2. • • ^ c] 

, and use an approximation ( V(T) ~ 

y(( T >) = y(T). Note that tiie equality ( T > = T is true because each 
of the components T,, for j = 1, C, is an independent random variable 
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whose uncertainty depends only on the probability distribution over the 
destinations of riding and waiting passengers assigned to car j. 

[037] For the same reason, this approximation is quite good on average. 
The exact landing time of each car T - depends, of course, on earlier 
assignments made to existing passengers, and their uncertain 
destinations. In other words, the landing pattem depends indirectly on the 
expected waiting time 21 1 of the existing passengers 111-113. A method 
for determining 210 the expected waiting time 21 1 of existing passengers 
1 1 1-1 14 is described by Nikovski et al., in '^Decision-theoretic group 
elevator scheduling ^ 13^ Intemational Conference on Automated 
Planning and Scheduling, June 2003, and U.S. Patent Application Sn. 
10/161,304 ''Method and System for Dynamic Programming of Elevators 
for Optimal Group Elevator Control filed by Brand et al. on June 3, 
2002, incorporated herein by reference. For short, this method is referred 
to as the "Empty the System Algorithm by Dynamic Programming" 
(ESA-DP) method. 

[038] So far we have considered the parking pattem T and T as 
functions of di fixed existing assignment of passengers to cars. However, a 
current decision of the scheduler 200, i.e., which car should be assigned 
to service the new passenger 113, changes this assignment. Because the 
scheduler can select any one of C cars, there are C possible resulting 
assignments, and hence, C possible distributions of the landing pattem 
219. If we use the above approximation, then we need the landing pattem 
^(0 = [^ib ^i2, .-M ^ic],for/= 1, C, 
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which occurs when the new passenger 113 is assigned to car i. The 
meaning of each entry T -j is the expected landing time of car j when the 
new passenger 1 13 is assigned to car i. 

[039] After the matrix for the landing pattem 219 for the C cars has been 
built, the expected cumulative waiting time 221 of future passengers 214 
corresponding to each of the landing pattem, i.e., rows of the matrix, can 
be determined. 



[040] We provide a procedure for determining an expected waiting time 
of future passengers 214 as a function of any landing pattem 219 
T = [Ti, T2, Tq]' 

[041] Because the waiting time 221 of the future passengers 214 is 
invariant with respect to the particular order of car arrivals, i.e., it makes 
no difference whether car "2" arrives in ten seconds and car "3" arrives in 
fifty seconds, or vice versa. We sort the landing pattem T 219 in an 
ascending order: 0 < Ti < 72 < . . . Tc. With this assumption, we define 
V^(T) as the expected cumulative waiting time 221 of all future 

passengers 114 within the time interval t ^ [0, Tq]: V^(T) = {^"KO^^, 

where n{t) is the expected number of passengers waiting at the main floor 
103 at time t. 



[042] Before describing our car assignment procedure, we introduce 
exponential discounting of future waiting times 221 because of a bias in 
the predicted parking times of the cars. The bias is due to our 
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approximating assumption that no future arrivals above the main floor 
occur before the end of the current landing pattem. 

[043] In practice, such future arrivals do occur, albeit infrequently. These 
passengers will be assigned to cars with riding and waiting passengers. 
Those cars are then delayed in reaching the main floor. Thus the landing 
times estimated by the ES A-DP process may underestimate slightly the 
actual times for near future predictions, and, perhaps, significantly for 
far-future predictions. 

[044] The near future can be defined as the average time it takes a car to 
make a round trip from the main floor and back, for example 40-60 
seconds for a medium sized building. This time is computable. 

[045] One way to discount estimates far into the future is to multiply the 
estimates by exp(-^0» where ^ > 0 is a discounting factor. 

[046] Similarly to the case above, we define the expected discounted 
cumulative waiting time of future passengers to be V^(T) = 



e'^'n{t)dt . The interval [0, Tq] can be split into C different intervals 

[T'i-b T\\, for / = 1, . . C, setting To = 0. The expected number of 
passengers waiting at time t ^ [Ti.i, T{\ is proportional to the time elapsed 
since the last time a car landed at the main floor was {Ti,\). 

[047] If we model the arrival of future passengers 1 14 as a Poisson 
process with a rate^ , then the expected number of passengers at the main 
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floor is n(t) = ^(t- Ti.{), and the integral above splits into C parts that can 
be evaluated. We assume that the cars can pick up all passengers waiting 
at the main floor instantaneously, because loading times are small 
relative to waiting times. 

[048] However, if car i reaches the main floor and finds it empty, then it 
does not depart immediately at its arrival time Ti. Instead, the car waits at 
the main floor until a future passenger 1 14 turns into the new passenger 
1 13 on signaling 120 a call. If there are j cars at the main floor at time t = 
0, then the first j passengers do not wait at all. Each passenger boards a 
car immediately, with no waiting time. The significant but speculative 
savings in this scenario are balanced against a real cost of not using those 
cars to service a new passenger at an upper floor. In order to quantify 
these savings, the elevator cars at the main floor are modeled accurately. 

[049] Semi-Markov IModel 

[050] To correctly estimate the waiting time 221 of future passengers 
214, given the actual behavior of cars when nobody is waiting at the main 
floor, we employ a semi-Markov chain whose states and transitions 
describe the behavior of cars landing at the main floor. 

[051] A semi-Markov chain includes a finite number of states 5,, i = 1, 
...,Ns, average momentary costs Uj , expected transition times ^ ,y , 
probabilities Py of the transitions between each pair of states 5,- and Sj , 
and an initial distribution ^(5,), which specifies the probability that the 
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system starts in state 5,, see Bertsekas, ''Dynamic Programming and 
Optimal Control'' Athena Scientific, Belmont, Massachusetts, 2000. 
Volumes 2, pages 261-264. Furthemiore, each semi-Markov chain 
contains an embedded fuUy-Markov chain evolving in discrete time, 
whose cumulative transition costs Rij are defined as Rij = ^ ifij , and all 
transitions are assumed to occur within a unit of time. The states in the 
semi-Markov chain used for our problem are labeled by the triple (ij, 
m), where i is the number of cars yet to land at the main floor, 7 is the 
number of cars parked currently at the main floor waiting for passengers, 
and m = C - / - 7 is the number of cars already departed from the main 
floor. 

[052] As shown in Figure 3, we organize the states of the semi-Markov 
chain in a two-dimensional grid or matrix. Each element Sim 301 in the 
matrix 300 corresponds to a state (ij, m). The grid structure in Figure 3 
is for an embedded semi-Markov chain for a building with four shafts. 
Row 302 / of the model contains all possible states of the system just 
after car i has arrived at time and has picked up all passengers that 
might have been waiting at the main floor. Note that the vertical time axis 
303 is not drawn to scale. Only transitions shown in bold arrows 304 
have non-zero costs. The cost of all other transitions is zero. Transitions 
labeled with n+ 305 for some number n are taken when n or more 
passengers arrive. 

[053] First, we provide a solution for the generic situation represented by 
this model, namely when no cars are parked at the main floor at the 
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current decision time (Ti > 0), and later extend the solution to the case 
when some cars are parked at the main floor. 

[054] For the generic case, the starting state of the chain is a state (C, 0, 
0), i.e., all C cars are yet to land at the main floor. The terminal states are 
those in the bottom row of the model, when all C cars have landed, and 
depending on how many of the future passengers have arrived in the 
interval t ^ [0, Tc\. Either all cars have departed with passengers on 
board, i.e., state (0, 0, C) 210, or some cars are still present at the main 
floor, i.e., states (0,7, C - j) for some j > 0. 

[055] Each state (ij, m) in the rows above the bottom one (i > 0), where 
j=C-i-m, can transition to two or more successor states. This depends 
exactly on how many future passengers arrive during a time interval t ^ 
[Ti, Ti+i]. For example, the chain transitions from state (4, 0, 0) to state 
(3, 1, 0) only when no passengers arrive by time Ti, and transitions to 
state (3, 0, 1) when one or more passengers arrive by that time. Each of 
the transitions in Figure 3 is labeled with the number of passengers that 
should arrive when this transition is taken. 

[056] The time to con^lete each transition is readily determined to be 
the interval ^Ti = Ti- Tf.i between the arrival of two cars. The 
probabiUty of each transition can also be determined because the 
transition is equal to the probability that a particular number of future 
passengers arrive within a fixed interval from a Poisson process with 
arrival rate . Thus, the probability p(x) that exactly x passengers arrive 
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in time A Ti is p(x) = (A&T.ye'^^^ lx\. For transitions labeled with an exact 
number of arriving passengers, this formula can be used directly. For 
transitions labeled with n+, meaning that they are taken when n or more 
new passengers arrive, the probability of the transition is one minus the 
sum of the probabilities of all remaining outgoing transitions from this 

state: p{n+) = 1 - S^-o^W- 

[057] Determining the cost of transitions labeled with an exact number of 
passengers is straightforward because the number of arriving passengers 
is less than or equal to the number of cars parked at the main floor. None 
of these passengers has to wait, and the cost of the corresponding 
transitions is zero. However, determining the cost of the last or rightmost 
transition from each state is quite involved. Such a transition corresponds 
to the case when n or more passengers arrive at the main floor, while 
only n - 1 cars are parked there. The computation has to account for the 
fact that if x future passengers arrive, and x-n, the first n - 1 of 
passengers take a car and depart without waiting, and only the remaining 
x-n+\ passengers have to wait. 

[058] Figure 3 shows that for any state Sim of the grid, as defined above 
and j=C -i -m, the transition shown in bold is taken when more than j 
future passengers arrive, i.e., n=j - 1. Hence, if that transition is taken 
and X future passengers arrive, then only the last x-j passengers have to 
wait. In other words, if x passengers appear within some time t, the 
differential or momentary cost at that time isx- j. 
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[059] Because such a transition covers the case when some number of 
passengers greater than j appear, and this number can theoretically be 
arbitrarily large, even in a finite time interval, the expected cost of the 
transition is a weighted sum over all possible numbers of arrivals x, from 
7 + 1 to infinity, and the weights are the probabilities that x arrivals occur, 
as given by the Poisson distribution. 



[060] In addition, the differential costs at time t can be discounted by a 
factor of exp(-^f)» as described above. This reasoning yields the 
following expression for the expected discounted cumulative waiting 
time im of main floor passengers during the last transition out of state 
Sim, with j = C'i-m: 



[061] After a change of integration variables, simplification, and splitting 
of the integral into two parts according to the two components of the 
difference between x - j, the expression for the cost evaluates to 
R/3.^ = e^ ^-'[F( A Tc,m) - F(0)], making use of a function 

for some arbitrary, but fixed integration constant Cq, which we set to zero 
for convenience. 
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[062] After all costs and probabilities of the semi-Markov model have 
been determined as described above, the cumulative cost of waiting 
incurred by the system when it starts in any of the model states can be 
determined efficiently by means of dynamic programming, starting from 
the bottom row of the model and working upwards, see Bertsekas, 
'''Dynamic Programming and Optimal Control^ Athena Scientific, 
Belmont, Massachusetts, 2000, Volumes 1, pages 18-24. Because the 
states in the bottom row are terminal and mark the end of the landing 
pattem, we set their waiting times to zero, i.e., we are not interested in 
the amount of waiting time accumulated after the last landing. 

[063] After the waiting times for all states are determined, we can obtain 
the cumulative waiting time for the entire pattem T from the initial state 
of the model. In the generic case, if there are no cars at the main floor at 
time ? = 0, then the initial state is always (C, 0, 0). The special case, when 
one or more cars are parked at the main floor at time ^ = 0, can be 
handled just as easily. In this special case, the starting state is (C - 1, 1, 0), 
where 1 is the number of cars at the main floor, and the expected 
discounted cumulative wait for the entire pattem is the waiting time of 
this starting state (5'c.i,o)- This eliminates the need to handle this special 
case separately from the generic one. 

[064] The procedure described above provides estimates V\P^vP(Ii) of 
the expected cumulative discounted waiting time 221 of future 
passengers 114, based on each of the landing pattem T, 219 resulting 
from the decision to assign the current call 201 to car i, i = 1 ..C. 
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Simultaneously, the ESA-DP process in step 210 determines estimates Wi 
of the cumulative non-discounted waiting time 21 1 of the existing 
passengers 21 1-213, including the new passenger 213 that signaled the 
call 201, when the call is assigned 230 to car i, / = 1, . . C. 

[065] In order to arrive at an optimal decision that balances the wait 211 
of existing passengers and the wait 221 of future passengers, the two sets 
of values and W/ are combined 230 to determine the adjusted waiting 
time 231. 

[066] There are significant differences between these two measures: The 
cumulative waiting time 21 1 of passengers W/, i.e. waiting 112 and the 
new passenger 213, is not discounted, while the cumulative waiting time 
221 of the future passengers 214 is discounted. 

[067] Furthermore, an objective of the scheduling process 200 is to 
minimize an average waiting time, and not the cumulative waiting time 
over some interval. For the purposes of optimization, the two measures 
are interchangeable only when the time intervals for all possible decisions 
are equal. 

[068] In general, this is not the case. The landing pattem does not have 
the same duration for each car. Therefore, the scheduling process 200 has 
to average waiting times from their cumulative counterparts. 
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[069] Obtaining the average expected waiting time 21 1 of existing 
passengers 11-113 from the cumulative waiting time is 

straightforward. The number N of existing passengers 1 1-1 13 is always 
known by the scheduler and does not depend on the candidate car number 
i, so ^W^/ N , On the other hand, obtaining the average waiting time 
221 of future passengers 214 from the cumulative discounted waiting 
time Vi^ over the duration of a landing pattem 219 is not as obvious. 

[070] The duration Tc of the landing pattem is known. If the arrival rate 
at the main floor is ^ , then the expected number of arrivals within Tc 
time units is^ Tc. However, dividing V, by ^ Tc is meaningless, because 
Vi has been discounted at a discount rate ^, 

[071] Instead, the discount factor exp(-^0 is an averaging weight for 
time t. If n(t) is the expected momentary number of passengers arriving at 
time t, as reflected in the costs of the Markov model, then Vj^ = 

e n{t)dt means the expected cumulative weighted number of 
passengers arriving during the time interval [0, Td- Therefore, the 
quantity ^ f e^n{t)dt/ ^ e^dt is the expected average number of 
future passengers arriving within this interval, properly normalized by the 
integral sum of all weight factors. Furthermore, Little's law specifies that 
w = ^ ^ /, see Cassandras et al., '"^Introduction to discrete event systems'^ 
Kluwer Academic PubUshers, Dordrecht, The Netherlands, 1999. This 

finally yields the time-normalized expected wait of future passengers 221 
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[072] Having obtained comparable estimates 211 and m 221 of the 
waiting times of existing and future passengers, these waiting times are 
combined 230 into a single adjusted waiting time 2231, for example by 

means of a weight 0-^-1, such that the adjusted waiting time is ^ ^ 

+ (l-«)^. 



[073] The balance between existing and future waits depends on how 
quickly the system can free itself of present constraints by delivering 
passengers. 

[074] Thus the optimal value of ^ can be determined empirically based 
on physical operating characteristics of the elevator system. We find that 
weight values in the interval [0.1, 0.3] stably produce acceptable results, 
regardless of the height of the building and number of shafts. 



[075] Effect of the Invention 



[076] The system and method as described herein can significantly 
reduce waiting time with respect to the conventional scheduling 
processes, with savings in the range of 5%-55%. These improvements are 
attributed to the look-ahead policy for future passengers. Elevator 
performance in up-peak traffic typically determines the number of shafts 
a building needs. Using standard guidelines for fitting elevators in a 
building, the invention can often reduce the number of required shafts for 
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mid- and high-rise office buildings by one, while still providing superior 
service. For a medium sized building, e.g., 25-30 floors, the cost per 
elevator can be about $200,000. Eliminating a shaft not only reduces the 
cost of the building but also the cost of maintenance, while increasing 
usable floor space. 

[077] Although the invention has been described by way of examples of 
preferred embodiments, it is to be understood that various other 
adaptations and modifi^cations may be made within the spirit and scope of 
the invention. Therefore, it is the object of the appended claims to cover 
all such variations and modifications as come within the true spirit and 
scope of the invention. 



23 



